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Abstract 

As part of the Wallops Coastal Oceans Observing Laboratory 
(Wa-COOL) Project, we sampled a time-series transect 
in the southern Mid-Atlantic Bight (MAB) biweekly. Our 
2-year time-series data included physical parameters, nutri- 
ent concentrations, and chlorophyll a concentrations. A 
detailed phytoplankton assemblage structure was examined 
in the second year. During the 2-year study, chlorophyll a 
concentration (and ocean color satellite imagery) indicated 
that phytoplankton blooms occurred in January/February 
during mixing conditions and in early autumn under strati- 
fied conditions. The chlorophyll a concentrations ranged 
from 0.25 pg f to 15.49 pg P 1 during the 2-year period. 
We were able to discriminate approximately 116 different 
species under phase contrast microscopy. Dominant phy- 
toplankton included Skeletonema costatum, Rhizosolenia 
spp., and Pseudo-nitzschia pungens. In an attempt to deter- 
mine phytoplankton species competition/succession within 
the assemblage, we calculated a Shannon Weaver diversity 
index for our diatom microscopy data. Diatom diversity was 
greatest during the winter and minimal during the spring. 
Diatom diversity was also greater at nearshore stations than 
at offshore stations. Individual genera appeared patchy, with 
surface and subsurface patches appearing abruptly and per- 
sisting for only 1-2 months at a time. The distribution of 
individual species differed significantly from bulk variables 
of the assemblage (chlorophyll a) and total phytoplankton 
assemblage (cells), which indicates that phytoplankton spe- 
cies may be limited in growth in ways that differ from those 
of the total assemblage. Our study demonstrated a highly 
diverse phytoplankton assemblage throughout the year, with 
opportunistic species dominating during spring and fall in 
response to seasonal changes in temperature and nutrients in 
the southern MAB. 

Keywords: coastal; diversity; Mid-Atlantic Bight; 
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Introduction 

An increasing amount of research has focused on the response 
of marine ecosystems to alterations in phytoplankton diversity 
resulting from global climate change (Moisan et al. 20 1 1 , 20 1 2). 
Little is known about how the phytoplankton assemblage, at the 
species level, responds to human-induced perturbations. There 
is a need for long-term time-series data from coastal areas where 
there is significant human impact to understand how climate 
affects marine ecosystems. Many large-scale programs, includ- 
ing the Joint Global Ocean Flux Study (JGOFS), Bermuda 
Atlantic Time-Series Study (BATS), and Hawaiian Ocean Time 
Series [HOTS; (Chavez et al. 1990, Malone et al. 1993, Karl 
et al. 2001)], have addressed spatial and temporal variability 
over long time periods and have emphasized that carbon flow 
is mediated by phytoplankton species composition and produc- 
tivity. The National Oceanic and Atmospheric Administration’s 
Integrated Ocean Observing System (NOAA’s IOOS) program 
vision was to create “a fully integrated ocean observing system” 
in the coastal regions to understand the relationship between cli- 
mate and ecosystems for “sustained living marine resources” 
(NOAA 2007). One of our goals was to determine the under- 
lying parameters that drive changes in phytoplankton assem- 
blage composition and, therefore, carbon flow to better predict 
the effects of climate change. Temperature, a critical variable 
that determines phytoplankton succession as well as light and 
nutrients, affects phytoplankton diversity and perturbations in 
assemblage structure (Moisan et al. 2002, Smayda et al. 2004). 

Our study focuses on phytoplankton dynamics near the 
Delmarva Peninsula, which is located in an area under the 
influence of Chesapeake and Delaware Bays on the east coast 
of the USA. Previous studies and satellite imagery do not sug- 
gest that our study area is significantly impacted by either of 
these bays. There have been excellent descriptive studies of 
phytoplankton assemblage composition in the coastal region 
of the Delmarva Peninsula (Marshall 1976, 1984, Marshall 
and Cohn 1 987). Studies in Chesapeake and Delaware Bays 
have extensively assessed nutrient dynamics and primary pro- 
ductivity in relation to hydrographic conditions, and spatial 
and temporal variability in phytoplankton assemblage struc- 
ture (see Watling et al. 1979, Harding et al. 1986, Marshall 
and Lacouture 1986, Fisher et al. 1992, Roman et al. 2005). 
The dynamics of the coastal phytoplankton assemblage are 
less known in relation to nutrient availability, vertical struc- 
ture, and physical influences. However, assessment of the 
Delmarva coastal region is important to: 1) develop a base- 
line understanding of the assemblage for long-term climate- 
related assessments, 2) improve productivity estimates, and 
3) improve our understanding of the influence of phytoplank- 
ton diversity on biomass and primary productivity. 
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The study area was in the southern Mid-Atlantic Bight 
(MAB), off the coast of Assateague Island (Virginia); it is 
located near the NASA Goddard Space Flight Center Wallops 
Flight Facility. The Gulf Stream flows northward further off- 
shore. Temperate species are dominant in the region, with 
subtropical species sometimes appearing in warm core rings 
originating from the Gulf Stream (Cox and Wiebe 1979). 
Waters in the region are typically well-mixed in the winter 
months, with strong stratification in the summer. 

We investigated physical, chemical, and biological param- 
eters at five stations along a transect extending from approxi- 
mately 9.6 km to 40.7 km offshore (Figure 1, Table 1). The 
stations ranged in depth from approximately 9 m to 32 m 
(Table 1). Phytoplankton species were examined at Stations 
1 and 5 to analyze community structure and to gain a bet- 
ter understanding of the forces affecting the phytoplankton 
assemblage. We also examined the vertical structure of the 
assemblage to assess the intrusion of high-nutrient waters, the 
influence of irradiance, and the occurrence of a bottom chlo- 
rophyll a maximum. Biological, physical, and nutrient data 
were obtained from a series of cross-shelf transect cruises 
from June 2005 to December 2007 as a part of the Wallops 
Coastal Ocean Observation Laboratory (Wa-COOL) Project. 
Cruises were conducted bi-weekly to monthly, weather per- 
mitting, aboard the Marine Science Consortium’s RV Phillip 
N. Parker. Stations ranged from approximately 9.6 km (Station 
1) to 40.7 km (Station 5) offshore from Assateague Island 
(Virginia) in the southern MAB (see Figure 1 and Table 1). 

Although we examined the physical parameters, nutrient 
patterns, and biological characteristics at all five stations, the 
phytoplankton assemblage structure data were acquired at 
only Stations 1 and 5. These stations were chosen because 
of the distance from each other and include coastal and oce- 
anic phytoplankton communities. At each station, the same 
discrete depths were sampled on each cruise (Table 1). In 
addition to comparing assemblages immediately offshore 
vs. those further offshore, sampling various depths allowed 



Longitude (°) 

Figure 1 Cruise transect stations (red) offshore from Assateague 
Island. Stations are numbered, with Station 1 closest to coastline, and 
proceed to Station 5 furthest offshore. 


Table 1 Station information. 


Station 

Latitude 

(N) 

Longitude 

(W) 

Distance 

from 

shore (km) 

Total 

depth 

(m) 

Depths 
sampled (m) 

1 

37.8164 

75.3753 

9.6 

9.0 

0, 4,7 

2 

37.7913 

75.3144 

15.7 

12.0 

0, 6,9 

3 

37.7697 

75.2449 

22.2 

21.0 

0, 5, 10, 15 

4 

37.7386 

75.1507 

31.3 

23.0 

0, 5, 10, 20 

5 

37.6991 

75.0543 

40.7 

32.0 

0, 5, 10, 20, 29 


assessment of the vertical components of phytoplankton 
assemblage structure. 

The phytoplankton assemblage was dominated by diatoms 
during the spring and fall blooms, with high diversity during 
the summer. Nutrients never became limiting in the area, but 
there were other competing sources of carbon during summer, 
as demonstrated by high Redfield ratios. Grazing appeared to 
control the taxonomic composition of the assemblage, espe- 
cially with respect to Cyanobacteria (Moisan et al. 2010, J. 
Ambler et al. unpublished observations). Despite strong strat- 
ification, we often observed a well-developed chlorophyll 
maximum above the sandy bottom where nutrients were low. 

Materials and methods 

Depth-specific water samples were collected either by bucket 
(surface samples) or Niskin bottle, and held at in situ temper- 
atures until processing (Table 2). Upon return to the labora- 
tory, 20 ml samples were filtered under negative pressure <5 
mm Hg, preserved with 2% glutaraldehyde (final concentra- 
tion), stained with proflavine (final concentration 5 j_ig ml' 1 ), 
and immediately filtered onto 2.0 pm black polycarbonate 
filters (GE Osmonics, Swedesboro, NJ, USA), using 8.0 pm 
nitrocellulose backing filters to ensure an even distribution of 
cells during filtration (modified from Sherr et al. 1 993, Kemp 
et al. 1993). Polycarbonate filters were placed on moistened 
microscope slides, a drop of immersion oil was added, and a 
coverslip was added. Prepared slides were stored in a freezer 
until analysis. 

Phytoplankton assemblage composition by 
epifluorescence microscopy 

Phytoplankton taxa from net tows (mesh size 20 pm) were 
preserved and used to identify taxa >20 pm, resolved to 

Table 2 Seasonal mean fluorometric chlorophyll concentrations 
(±SD). 


Winter Spring Summer Autumn 

Station 1 6.40 (±1.44) 4.23 (+2.68) 5.51 (±3.86) 5.81 (±1.71) 

Station 2 4.88 (±1.61) 3.16 (±2.39) 3.21 (±3.33) 4.22 (±1.22) 

Station 3 4.61 (±2.88) 2.21 (±1.24) 2.07 (±1.53) 3.72 (±2.65) 

Station 4 3.89 (±1.86) 1.98 (±1.11) 1.49 (±1.30) 1.80 (±1.08) 

Station 5 3.26 (±1.93) 1.93 (±1.30) 1.43 (±1.65) 1.95 (±1.57) 

Dimensions are pg l" 1 (n=3). 
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genus or species level. Microscopic analysis was performed 
under an Olympus (Tokyo, Japan) BX-5 1 compound micro- 
scope with a mercury light source and a fluorescence filter 
cube (Olympus 41015) for proflavine fluorescence. An ocular 
grid was used to count ten fields on each slide, and fields were 
systematically chosen to cover all areas of the filter to take 
into account potentially uneven distributions of cells. Fields 
were identified, and initial counts and identifications were 
performed under 400x magnification. For small cells within 
the field, an oil immersion lOOx objective lens was used as 
well. When very few cells were included in ten fields, addi- 
tional fields were counted. All quantitative counts were made 
with whole water samples and were counted on filters using 
epifluorescence microscopy. Net tows were used only as a 
qualitative measure. We utilized different methods to reduce 
our error due to uneven distribution of cells and found that 
our % error results were within a coefficient of variation of 
20% (Moisan et al. 2010). 

Individuals that could not be identified down to the genus 
level were assigned unique four-digit identification numbers 
for classification purposes following Fryxell and Hasle’s con- 
ventions (Tomas 1997). Unknown diatoms were also classified 
as either centric or pennate. Known and unknown dinoflagel- 
lates were assigned to distinct size categories (<20 |im, 20-40 
pm, and >40 pm). With the exception of a few obvious spe- 
cies, flagellates were pooled together in a “flagellates” cate- 
gory because we did not have access to electron microscopy, 
which is needed for identifying these very small organisms. 
In general, we found high diversity, with approximately 116 
species that we were able to resolve with our methodology 
(Table 3). 

Conductivity, temperature, depth (CTD) 

A Sea-Bird 49 CTD (Bellevue, WA, USA) integrated into 
an optical profiling package was deployed at all stations on 
transect cruises to measure conductivity and temperature 
through all depths of the water column. For data collection, 
the profiler was lowered to the water surface for at least 2 
min to allow ample time for equilibration. The profiler was 
then lowered to 2 m above the seafloor and raised at a speed 
of approximately 0.5 m s 1 while logging data at a rate of 1 
Hz. Once back in the laboratory, CTD data were extracted 
from the profiler’s WetLabs (Philomath, OR, USA) DH-4 
data handler and processed. We converted from pressure 
(decibar) to depth (meters) following Sea-Bird’s Application 
Note 69 (from UNESCO 1983). Density was calculated 
using the UNESCO 1983 (EOS 80) polynomial. Finally, all 
data were binned at 1-m intervals and used to create water 
column profiles. 

Fluorometric chlorophyll a 

Phytoplankton counts were also compared with chlorophyll a 
data collected on each transect cruise. Water samples for fluo- 
rometric chlorophyll a and microscopic measurements were 
taken from the same Niskin bottles to ensure that we were sam- 
pling a consistent phytoplankton assemblage. Samples were 


collected on 0.7 pm Whatman GF/F filters (Piscataway, NJ, 
USA), stored in Histoprep tissue capsules (Fisher Scientific, 
Waltham, MA, USA), and flash frozen in liquid nitrogen until 
processing. For processing, each sample was extracted in a 
test tube with 7 ml of 90% acetone for approximately 21-24 h 
in a dark freezer. Chlorophyll a fluorescence was then mea- 
sured using a Turner Model 10-AU fluorometer (Sunnyvale, 
CA, USA) according to the method of Welschmeyer (1994). 

Dissolved inorganic nutrients 

Samples for nutrient analyses were collected by filtering water 
through 0.2 pm capsule filters (Pall 12140, (Port Washington, 
NY, USA) using a peristaltic pump. The filtrate was collected 
and immediately frozen until analysis at the University of 
Maryland Center for Environmental Science. The nutrients 
measured and their associated detection limits (in parenthe- 
ses) included ammonium (0.21 pM), nitrite (0.03 pM), nitrate 
+nitrite (0.01 pM), phosphate (0.02 pM), and silicate (0.18 
pM). Analyses were performed using a modified version of the 
method in the Chesapeake Biological Laboratory’s Technical 
Report Series No. 158-97 (Keefe et al. 2004). 

Particulate organic carbon and nitrogen (POC 
and PON) 

POC and PON levels were measured through August 2007. 
Whatman GF/F filters were pre-combusted in a muffle fur- 
nace at 450°C for 4-6 h prior to collecting samples. Water 
samples were then filtered (<5 mm Hg negative pressure) 
onto the GF/F filters, and frozen until analysis. Prior to anal- 
ysis, filters were dried at 50°C overnight in the laboratory. 
Once dried, they were taken to Horn Point Laboratory for 
final processing. The filters were put in a desiccator, fumed 
with HC1 to remove the inorganic carbon, and re -dried. 
Samples were then analyzed with an Exeter Analytical, Inc. 
CE-440 Elemental Analyzer (Buckhurst Hill, UK) following 
methods in the Horn Point Laboratory Standard Operating 
Procedures. 

Results 

We sampled every 2 weeks at five stations normal to the 
coast and gently sloped from 7 m to 32 m (F igure 1 ) . Shallow 
stations were well mixed to the bottom throughout the year. 
Stations 4 and 5 became stratified during the period from 
July to October, when primary productivity was highest. 
We conducted detailed phytoplankton microscopic analy- 
ses at the shallowest station (Station 1) and offshore Station 
5 for 1 year. During the study, measurements of physical 
and chemical characteristics of these sites were made for 
later interpretations of coincident biomass and biodiversity 
observations. 

Physical parameters 

Our two shallow inshore stations were generally well 
mixed to the bottom throughout the year and showed little 
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Table 3 Phytoplankton identified at Stations 1 and 5 from November 2006 to December 2007. 


Diatoms 

Achnanthes sp. 

Actinoptychus senarius (Ehrenberg) 

Asterionella glacialis (Castracane) Round 
Bacillaria paxillifera (Mueller) Hendey 
Bacteriastrum deliculatum Cleve 
Biddulphia alternans (Bailey) Van Heurck 
Cerataulinci pelagic (Cleve) Hendey 
Chaetoceros sp. 

Chaetoceros aequatorialis Cleve 
Chaetoceros affinis Launder 
Chaetoceros compressus Lauder 
Chaetoceros costatus Pavillard 
Chaetoceros curvisetus Cleve 
Chaetoceros danicus Cleve 
Chaetoceros debilis Cleve 
Chaetoceros decipiens Cleve 
Chaetoceros diadema (Ehrenberg) Gran 
Chaetoceros didymus Ehrenberg 
Chaetoceros lorenzianus Grunow 
Chaetoceros socialis Lauder 
Corethron criophilum Castracane 
Coscinodiscus radiates Ehrenberg 
Coscinodiscus wailesii Gran et Angst 
Cyclotella litoralis Lange et Syvertsen 
Dactyliosolen fragilissimus (Bergon) Hasle 
Detonula sp. 

Ditylum brightwellii (West) Grunow 
Eucampia zodiacus Ehrenberg 
Fragilaria sp. 

Gonioceros sp. H. et M. Peragallo 
Guinardia flaccida (Castracane) H. Peragallo 
Guinardia striata (Stolterforth) Hasle 
Heliotheca sp. Yankovskii 
Dinoflagellates 

Dinoflagellates <20 pm 
Gonvaulax sp. Diesing 
Gymnodinium sp. Stein 
Gyrodinium sp. Kofoid et Swezy 
Katodinium rotundatum (Lohmann) Fott 
Amphidinium sp. Claparede et Lachmann 
Unknown #3014 
Unknown #3023 
Dinoflagellates 20-40 pm 
Amphidinium sp. Claparede et Lachmann 
Amphidinium sphenoides Wiilff 
Gonyaulax polygramma Stein 

Gonvaulax spinifera (Claparede et Lachmann) Diesing 
Gymnodinium sp. Stein 
Gyrodinium sp. Kofoid et Swezy 
Prorocentrum gracile Schiitt 
Protoperidinium divergens Ehrenberg 
Dinoflagellates >40 pm 
Alexandrium sp. Halim 
Ceratium fusus (Ehrenberg) Dujardin 
Ceratium lineatum (Ehrenberg) Cleve 
Ceratium macroceros (Ehrenberg) Vanhoffen 
Ceratium trichoceros (Ehrenberg) Kofoid 
Ceratium tripos (Muller) Nitzsch 
Dinophysis acuminata Claparede et Lachmann 
Dinophysis dens Pavillard 
Dinophysis norvegica Claparede et Lachmann 
Gymnodinium sp. 

Prorocentrum micans Ehrenberg 


Hemiaulus sinensis Greville 

Leptocylindricus minimus Gran 

Lioloma c.f. pacificum (Cupp) Hasle 

Melosira sp. Agardh 

Navicula transitans Cleve 

Nitzschia longissima (Brebisson) Ralfs 

Odontella mobilensis (Bailey) Grunow 

Odontella sinensis (Greville) Grunow 

Paralia sulcata (Ehrenberg) Cleve 

Plagiogramma sp. Greville 

Pleurosigma normanii Ralfs 

Pseudo-nitzschia multiseries (Hasle) Hasle 

Pseudo-nitzschia pseudodelicatissima (Hasle) Hasle 

Pseudo-nitzschia pungens (Grunow ex Cleve) Hasle 

Rhabdonema sp. Kiitz 

Rhaphoneis amphiceros Ehrenberg 

Rhizosolenia sp. (Ehrenberg) Brightwell 

Rhizosolenia alata Brightwell 

Rhizosolenia imbricate Brightwell 

Rhizosolenia setigera Brightwell 

Rhizosolenia styliformis Brightwell 

Skeletonema costatum (Greville) Cleve 

Thalassionema nitzschioides (Grunow) Mereschkowsky 

Thalassiosira sp. Cleve 

Thalassiosira (20 pm) 

Thalassiosira (30 pm) 

Thalassiosira (40 pm) 

Thalassiosira (50 pm) 

Thalassiothrix sp. 

Unknown pennate #1024 
Unknown pennate #1037 
Unknown centric #1040 
Unknown pennate #1061 


Unknown #3027 
Unknown #3033 
Unknown #3035 
Unknown #3036 
Unknown #3039 
Unknown #3045 


Protoperidinium pellucidum (Bergh) Schiitt 

Pyrophacus horologicum Stein 

Scripsiella trochoidea (Stein) Loeblich 

Unknown #3015 

Unknown #3016 

Unknown #3018 

Unknown #3021 

Unknown #3053 

Prorocentrum scutellum Schroder 
Protoperidinium sp. 

Protoperidinium depressum (Bailey) Balech 

Protoperidinium pallidum (Ostenfeld) Balech 

Unknown #3040 

Unknown #3042 

Unknown #3044 

Unknown #3046 

Unknown #3051 

Unknown #3056 


See Table 1 for station locations. 
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variability, whereas the outer shelf stations generally strati- 
fied during the summer and early autumn (Figure 2). Water 
temperature in the region ranged from ~5°C in March to 
~26°C in September during the 2-year period. Salinity 
ranged from 30.76 to 33.55 during the study period, and 
we observed a gradient of higher salinity with distance 
offshore. 

Irradiance 

Photosynthetically available radiation reached the bottom at 
the inshore stations (Stations 1 and 2), and phytoplankton 
generally did not develop a depth-related chlorophyll maxi- 
mum because the mixed layers reached down to the bottom 
(data not shown). Station 5 appeared to develop a chloro- 
phyll maximum that was associated with low light levels, low 
micronutrients, and punctuated periods of stratification, espe- 
cially during spring and fall blooms. 


Fluorometric chlorophyll a 

Fluorometric chlorophyll a concentrations ranged from 0.2 
pg l' 1 to 15.5 pg F 1 over the course of the time series. Values 
were generally highest at Station 1 and declined with distance 
from shore (Figure 3). Mean surface chlorophyll a concentra- 
tions ranged from 1.6 (±1.3) pg F' at Station 5 to 4.6 (±3.0) 
pg F 1 at Station 1 . When examining the time series over the 
2-year period, it was evident that spikes in chlorophyll a 
were more common inshore throughout the year. It was also 
common to observe chlorophyll a maxima near the bottom, 
especially at Station 5, where the bottom chlorophyll a con- 
centration was over 29 times higher than the surface concen- 
tration in September 2007. 

During the spring bloom, we repeatedly observed a 
stratified pattern relating to temperature and salinity. 
Chlorophyll a values were higher at near-bottom depths. 
Nutrient concentrations were relatively high throughout 
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Figure 2 Temperature and salinity time-series data for Stations 1,3, and 5. Values are from continuous conductivity, temperature and depth 
(CTD) samples and binned at 0.5 m intervals. Values represent the median between the upcast and the downcast of the CTD. See Table 1 for 
station information. 
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Figure 3 Average fluorometric chlorophyll a concentrations and carbon to chlorophyll a ratio (C:Chl a) time-series plots for Stations 1, 3, 
and 5 (n=3). See Table 1 for station information. 


the water column during this period of time. It appeared 
that Station 1 was able to maintain a critical biomass for 
a longer period of time, probably due to less light limita- 
tion. Highest average concentrations of chlorophyll a typi- 
cally occurred in the winter, followed by the autumn bloom 
(Table 2). 

In contrast, we observed short bursts of certain taxa that 
may have been due to changes in water mass, upwelling/ 
downwelling events, or grazing. The increases in phyto- 
plankton biomass occurred in response to favorable upwell- 
ing conditions that occurred in small-scale patches along the 
coast (Wong 2002). It was difficult to determine whether 
these conditions were favorable because of nutrient concen- 
trations or light-availability conditions. Nevertheless, we 
postulate that a mixture of r- and k-selected types of phy- 
toplankton took advantage of these vertical mixing events 
(Silva et al. 2009). 


Particulate organic carbon and nitrogen 
(POC and PON) 

Throughout the study, POC ranged from 54 pg l' 1 (Station 
5) to 1563 pg T 1 (Station 1), and PON ranged from 10 pg 
f 1 (Station 5) to 221 pg T 1 (Station 1). Both POC and PON 
values typically decreased with distance from shore. POC 
values for Station 1 were highest during summer months at 
bottom depths but peaked in November 2006 at all depths. 
In November 2006, a concentration of 1563 pg l' 1 POC was 
measured at 7 m depth. At this inshore station, POC was low- 
est at the surface, and did not track any seasonal pattern. No 
consistent seasonal or vertical patterns for POC were evident 
at Station 5, but values ranged through an order of magnitude 
from 54 pg l' 1 to 534 pg l 1 . PON was typically highest at bot- 
tom depths at Station 1, and ranged from 47-221 pg l 4 . At 
the offshore station, PON ranged from 10 to 88 pg T 1 , but no 
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vertical pattern was apparent. No seasonal patterns for PON 
were observed at any station. 

The particulate C:N ratio ranged from 4.13 to 18.38. The 
C:N ratio of particulates was less than the Redfield Ratio of 
106:16 in 58% of samples (n=277). The particulate organic 
C:N ratio was less than the Redfield Ratio of 106:16 in 44% 
of the samples from Station 1 (n=57) and in 67% of samples 
from Station 5 (n=83). This may indicate that there was a sig- 
nificant amount of detritus in the area, which is considered 
labile. In addition, we observed a significant portion of detri- 
tal absorption (T.A. Moisan, data not shown). POC values 
were used to calculate the carbon to chlorophyll a (C:chl a) 
mass ratio. Patterns in the C:chl a ratio were similar across 
stations, but the values were more variable at stations further 
offshore, where short-duration, high intensity peaks occurred 
(Figure 3). During the late 2007 peak at Station 5, the C:chl a 
ratio was over 500. No clear seasonal, annual, or depth trends 
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in C:chl a ratios were apparent at any station. However, a 
multiple linear regression indicated that the C:chl a was sig- 
nificantly (R 2 =0.86) associated with detritus and temperature, 
based on detrital absorption (p<0.001) and temperature data 
(p=0.005). 

Dissolved inorganic nutrients 

Throughout the study, we measured micronutrients, including 
nitrate, phosphate, silicate, and ammonium, at select stations. 
Individual nutrients were highest in concentration at inshore 
stations and episodes of elevated concentration were consis- 
tently observed, as seen in the examples shown in Figures 4 
and 5. Nitrate, phosphate, and silicate had punctuated periods 
of elevated concentration at inshore stations, which dissipated 
with distance from shore. Nitrate peaked throughout the water 
column in December 2006 and November 2007. The mean 
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Figure 4 Micromolar concentrations of nitrate (left) and phosphate (right) at Stations 1, 3, and 5 (n=3). See Table 1 for station information. 
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Figure 5 Molar ratios of dissolved nutrients at Stations 1, 3, and 5 (n=3). See Table 1 for station information. 


nitrate concentration in December 2006 in the water col- 
umn at Station 1 was 2.21 pM (±0.078, n=3), and decreased 
offshore to a low of 0.46 pM (±0.261, n5) at Station 5. In 
November 2007, the mean nitrate concentration reached 3.07 
pM at Station 1 (±0.261, n=3). Elevated spikes of phosphate 
and silicate were observed during July through November 
2006 and were most evident at Stations 1 and 2. A second, 
more intense episode of phosphate and silicate elevation was 
noted in the autumn of 2007 at Station 1. The phosphate and 
silicate spikes was also noted at other stations, but occurred 
later in the season. 

Molar ratios of nitrogen, phosphate, and silicate indicated 
that this region is often nitrogen-limited (Figure 5). The 
nitrogen to phosphate molar ratio (N:P) typically remained 
less than the Redfield ratio (C:Si:N:P= 106: 15:16:1) at all 
stations. Nitrogen was a limiting nutrient in 97% (n=62) 
and 95% (n=86) of samples at Stations 1 and 5, respectively. 
Silicate also limited the growth of diatoms, as the Si:P ratio 


was less than 15:1 in 99% of samples (n=89); the silicate to 
phosphate (Si:P) ratio was highest at Station 1 and decreased 
offshore. 

Phytoplankton taxonomic composition 

Diatoms made up the majority of the phytoplankton biomass 
based on chlorophyll a, and typically accounted for about 40- 
50% of the total phytoplankton community. Phytoplankton 
concentrations ranged from 118 cells ml' 1 (November) to 
6104 cells ml' 1 (July) at Station 1. At Station 5, the range 
was roughly <100 cells ml' 1 (October) to 10,856 cells ml' 1 
(September). The phytoplankton community consisted of a 
range of taxa from small flagellates (~2. 5-5.0 pm) to large 
diatoms (>200 pm) to colonial Phaeocvstis. Together, the dia- 
toms Skeletonema costatum, Chaetoceros, Pseudo-nitzschia, 
Thalassiosira, Asterionella glacialis, and Rhizosolenia domi- 
nated throughout the study. 
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In contrast, dinoflagellates occurred at lower abundances 
and were most abundant in summer when waters were stratified 
(Figure 6). Dinoflagellates were common at the surface in late 
July and late August, with over 700 and 900 cells ml 1 , respec- 
tively. Dinoflagellates were most abundant at Station 5 during 
mid- June and ranged from 16 cells ml' 1 at the surface to 848 
cells ml' 1 (Figure 6). Throughout the year, ambient dinoflagel- 
late concentrations were generally <100 cells ml' 1 . Phaeocvstis 
colonies were present in the mid- to late summer. Cyanobacteria 
were at their highest concentrations in summer months and were 
associated with warm water (Moisan et al. 2010). 

Temporal distribution of the phytoplankton 
assemblage 

Individual phytoplankton taxon distribution was patchy both 
spatially and temporally. It was common to see a short burst 
of one species that would almost disappear from the assem- 
blage by the next cruise (~2 weeks later). There was tremen- 
dous variability in assemblage compositions by depth that 


resulted in high phytoplankton diversity. We saw short peaks 
of high concentration for certain taxa that may have been due 
to changes in water mass, up welling/ down welling events, or 
grazing (J. Ambler, unpublished data). We observed favor- 
able upwelling conditions, which occurred along the coast in 
small-scale patches (Wong 2002). 

Patterns of diatom diversity 

Throughout our study, we noted minor spatial and temporal 
variations in the diversity of the diatoms when comparing 
the assemblages at the inshore and offshore stations. Diatoms 
were classified down to the genus level and Shannon-Weaver 
diversity indices (FF) were calculated for each sampling date 
(Shannon and Weaver 1949, Pielou 1966). Diatom diversity 
ranged from 0.343 to 1.487 at the inshore station and from 
0.293 to 0.906 at the offshore station. Average diversity was 
lowest in the spring and autumn, when blooms typically 
occurred (Table 4). We chose several numerically abundant 
phytoplankton genera and present their seasonal distribution 
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Figure 6 Patterns of diatom, dinoflagellate, and flagellate cell concentrations at Stations 1 and 5. Note that these microscopic analyses began 
in November 2006. Coefficient of variation for each sample was <20%. See Table 1 for station information. 
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Table 4 Seasonal mean diatom genus diversity (Shannon Weaver index) (±SD). 


W inter 


Spring 


Summer 


Autumn 


Station 1 1.203 (±0.402) 0.568 (+0.110) 0.804 (±0.149) 0.753 (±0.206) 

Station 5 0.819 (±0.028) 0.538 (±0.107) 0.671 (±0.237) 0.595 (±0.191) 


here. The chosen genera had relatively high abundances and 
had unique and episodic temporal patterns throughout the 
year. Data for Thalassiosira, Pseudo-nitzschia, Rhizosolenia, 
are depicted in Figure 7. The appearances of taxa were mixed 
throughout the water and appeared to persist for more than 2 
weeks for each taxon. We also observed temporal patchiness 
in other genera, such as Chaetoceros, Skeletonema, and colo- 
nial Phaeocystis (Table 5). 

Discussion 

There has been an increasing emphasis on assem- 
bling long-term time-series data, including Long-Term 
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Ecosystem Research (LTER), CalCoFI (California 
Cooperative Oceanic Fisheries Investigations), BATS, 
HOTs, MOBY (Marine Optical Buoy), and NASA satel- 
lite products because they are ideally suited for investiga- 
tions that detect subtle habitat changes, irregularly spaced, 
stochastic events, and complex, interdependent ecological 
phenomena that affect biogeochemical cycles in the world 
ocean (Karl et al. 2001). The net flux of inorganic nutrients 
and organic carbon through the coastal zone is controlled 
by diverse biogeochemical processes that occur within a 
coupled web of ecological subsystems. We focus in this 
paper on a 2-year time series of chlorophyll a, physical 
parameters, and nutrients, with 1 year of detailed micros- 
copy of taxon composition. 
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Figure 7 Spatial distribution of Thalassiosira, Pseudo-nitzschia, and Rhizosolenia. Seasonal trends in cell concentrations at Station 1 (left) 
and Station 5 (right). Coefficient of variation for each sample was <20%. See Table 1 for station information. 
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Table 5 Abundances of dominant genera from November 2006 through December 2007 and peak depths. 


Genus 

Station 

Peak months 

Peak 
depth (m) 

Cell count at 
Peak depth 
(cells ml' 1 ) 

Mean cell count, 
all depths sampled 
(cells ml' 1 ) 

Thalassiosira 

1 

December (2006) 

4 

883 

655 (+276) 


5 

November (2006) 

5 

1580 

355 (+686) 

Pseudo-nitzschia 

1 

March 

4 

873 

691 (±225) 



August 

7 

608 

456 (±230) 



September 

7 

696 

448 (±337) 


5 

March 

5 

1285 

756 (±407) 



September 

29 

656 

307 (±249) 

Rhizosolenia 

1 

March 

0 

348 

255 (±81) 



September 

7 

920 

715 (±197) 


5 

March 

5 

275 

204 (±55) 



September 

5 

216 

96 (±76) 

Chaetoceros 

1 

August 

4 

1400 

928 (±586) 


5 

July 

10 

960 

435 (±401) 



September 

5 

896 

555 (±244) 



October 

20 

1032 

454 (±381) 

Skeletonema 

1 

December (2006) 

7 

787 

569 (±201) 



September 

4 

2480 

1803 (±628) 


5 

September 

29 

6136 

1578 (±2555) 

Phaeocystis 

1 

July 

4 

5112 

2403 (±2570) 


5 

November (2006) 

0 

6967 

3025 (±2571) 



July 

20 

1848 

566 (±834) 



September 

29 

2856 

571 (±1277) 


Cell concentrations had a coefficient of variation of <20% where n=19 depths. See Table 1 for station information. 


A limited number of studies have characterized the 
phytoplankton assemblage in Virginian coastal waters 
(Marshall 1991, Malone et al. 1991). The Virginian coastal 
ocean phytoplankton assemblage is considered diverse, 
often dominated by small, chain-forming diatoms and a 
picoplankton component of Cyanobacteria and chloro- 
phyceans (Malone et al. 1991). The annual cycle of phy- 
toplankton includes a peak in the spring, which is usually 
dominated by diatoms. By late spring, phytoplankters 
become nutrient-limited by either nitrate or phosphate and 
flagellates often dominate (Fisher et al. 1999). A bloom 
smaller in magnitude often occurs in the autumn and is 
different in composition from the spring bloom (Marshall 
1991). We found a diverse assemblage consisting of 116 
phytoplankton species dominated by different diatoms dur- 
ing the year, with pulses of flagellates during summer in 
response to localized increases in Cyanobacteria (Moisan 
et al. 2010). 

The phytoplankton assemblage was diverse and phyto- 
plankton succession was probably driven by a replacement 
of species by other species, as hypothesized by theories 
contributed by Smayda (1997, 1998). For example, we 
observed many patches of phytoplankton species over time 
that we believe were composed of opportunists that are able 
to take advantage of high nutrients in a mixed high-light 
environment. Assemblages were dominated by flagellates 
when nutrients were limiting in summer; heterotrophic 
flagellates are able to feed on Cyanobacteria (Moisan et al. 
2010). 


Seasonal patterns in phytoplankton biomass and 
taxonomic composition 

We were able to monitor physical properties, as well as chlo- 
rophyll a and POC consistently over the course of 2 years 
(Figures 2 and 3). In 2006, we found a weak winter/spring 
bloom, which was also confirmed by Ocean Color satellite 
imagery. In 2007, we were unable to capture the winter/spring 
bloom because of limited sampling opportunities, but satel- 
lite imagery indicated a stronger bloom, peaking in January 
and February. During both years (2006-2007), we observed 
a secondary bloom in the autumn, beginning in September 
and October in 2006 and 2007, respectively. Phytoplankton 
appeared to increase production in quick pulses in response to 
pulses in nitrate and phosphate, as shown by a great degree of 
variability in C:chl a ratios (Figure 4). High C:chl a ratios are 
probably representative of an association of detrital biomass 
and phytoplankton, as shown by the relatively high detrital 
matter to phytoplankton particle ratio (N. Ohi et al., unpub- 
lished data). It is likely that a portion of POC was contributed 
by detritus; we observed high detrital absorption in compari- 
son to phytoplankton absorption (see Results). 

It appears that the phytoplankton assemblage was able to 
acclimate to deep mixed layers because it was associated with 
deep euphotic zones that often reached the bottom at both sta- 
tions, with a range of chlorophyll-specific absorption values 
[a*ph (A,) associated with chromatic adaptation (T.A. Moisan, 
unpublished data)]. There was little relation between chloro- 
phyll a and nutrient ratios (Figure 5); the limiting nutrient in 
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this region appeared to be nitrate. Based on nutrient profiles, 
we postulate at least two sources for nutrients. The higher 
concentrations and increased chlorophyll a concentration at 
near-bottom depths of Station 5 led us to believe that there 
was a sediment-based flux of nutrients at this offshore sta- 
tion. However, the major nitrate spike in December 2006 and 
phosphate spike in July 2006, which were largest at Station 1 
and decreased with distance from shore, indicate that a sec- 
ond source of nutrients may be terrestrially-based. Because 
Chincoteague Bay has such a low flushing rate, it is not 
likely that it is the source of these nutrients during periods of 
punctuated growth. N. Ohi (unpublished data) found a tight 
correlation between phosphate and optical properties during 
the time series described by inherent optical properties and 
bulk reflectance. It is likely that nutrients were replenished 
by longshore flow from large sources of estuarine water from 
Delaware Bay and possibly the Chesapeake Bay or upwelling 
events (A. Kumar, pers. communication). 

Diversity patterns of the phytoplankton assemblage 
in the coastal environment 

To quantify phytoplankton community diversity, we utilized 
the Shannon Weaver index, which is a commonly used metric 
in ecology. Because we relied heavily on preserved samples 
enumerated using contrast microscopy, in addition to epi- 
fluorescence, our analysis of size and taxonomic diversity is 
constrained to those forms that were well preserved and iden- 
tifiable by light microscopy. Therefore, our Shannon Weaver 
index is most likely an underestimate of the diversity in our 
study area. It appeared that diversity indices were variable 
over the course of the year, with a decrease in diversity dur- 
ing the autumn bloom and greater diversity in the late winter/ 
early spring bloom. Additionally, we observed episodic peaks 
in certain taxa, which appeared to be associated with peaks in 
nutrients and possibly different water masses. 

We present six different taxa, including Thalassiosira, 
Pseudo-nitzschia, Rhizosolenia, that exhibited a patchy tem- 
poral distribution (Figure 7). A strong cohesiveness with 
temperature and the degree of stratification and turbulence 
were observed at Stations 1 and 5. Stratification and the depth 
to which water was generally mixed downward appeared 
to affect phytoplankton diversity but it was also indirectly 
affected by zooplankton diversity via changes in phytoplank- 
ton size distribution (J. Ambler, unpublished data). Model 
simulations in multi-species simulations suggest that global 
patterns of phytoplankton diversity are determined by a few 
dominating species within each group (Adjou et al. 2012). 
We observed short successional cycles, which were possi- 
bly dependent on coastal upwelling events, and patches of 
certain taxa appeared to persist from roughly 2 to 6 weeks 
(Figure 7). We postulate that the persistence of these patches 
was related to the duration and strength of upwelling events, 
nutrient availability, and punctuated changes in temperature 
(Spatharis et al. 2011). Dinoflagellates and smaller flagellates 
appeared to be relatively more abundant in late summer when 
nutrients had become depleted. It appears that the distribution 
of phytoplankton classes occurred along a gradient of a “r 


vs. K” growth strategy, where diatoms (r-selected) exploited 
well-mixed, turbulent, nutrient-rich conditions, whereas dino- 
flagellates (K-selected) dominated in stable, stratified waters 
with low nutrient concentrations (Silva et al. 2009). 

In summary, this study provides a 2-year time series that 
demonstrates how the MAB phytoplankton coastal assemblage 
responds to variability in wind forcing, which governs rela- 
tively short duration upwelling and downwelling events and 
temperature. The phytoplankton community was dominated 
by the diatoms Skeletonema costatum, Rhizosolenia spp., and 
Pseudo-nitzschia pungens throughout the year. Within the 
seasonal pattern, a series of short-term successions of various 
species, from r-strategists towards K-strategists, was associ- 
ated with hydrodynamic changes promoted by upwelling/ 
relaxation events (Silva et al. 2009). Our study detected high 
amounts of inter-annual variability in phytoplankton biomass 
and seasonal changes in phytoplankton diversity in the south- 
ern MAB. The data collected during this project will help to 
interpret marine biological changes and begin to distinguish 
between anthropogenic, climatologically forced, and natural 
plankton variability. 
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